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Abstract 

We propose a theoretical expression for the k- and ^-dependent dielectric function 
of a stack of two-dimensional layers coupled along the direction perpendicular to 
the layers, and discuss some of its properties. We argue that the plasma frequencies 
W) at k = should correspond to those which are experimentally obtained from optical 

measurements on e.g. La^^Sr^CuC^ via the /-sum rule analysis, regardless of the 
fact that such systems are strongly correlated. We discuss some of the ramifications 
due to strong anisotropy of the charge transport in these systems, and the lack of 
coherence for the transport in the direction perpendicular to the layers. 

oo 

In theoretical discussions of the high-T c cuprate superconductors it is usu- 
ally emphasized, that there exists a strong anisotropy in the electrodynami- 
cal properties between the in-plane and out-of-plane directions. Usually this 
^ anisotropy is regarded as resulting from a large mass anisotropy of the charge 

carriers, so that effectively m c ^> m a b where m c and m a b are the in-plane 
^ and out-of-plane effective masses. Ideally one prefers to use a model dielectric 

g function where these superconductors are treated as a stack of layers which are 

completely decoupled in the c direction, which amounts to taking m c — > oo. 
Such models are indicated as layered electron gas (LEG) models. An important 
consequence m c — > oo is, that the plasmon dispersion in the long- wavelength 
limit (kd <C 1 where k is the wavevector and d is the interlayer spacing) is 
given by u(k) ju^ k^j ' \jk 2 ah + k 2 c . For a given hnite value of k c this is an 
acoustical dispersion as a function of k a f,. If we write the dispersion in the form 
iL>(k)/uj a b sin#, where 6 is the angle between k and the c axis, it becomes 
clear, that k — > is a singular point in the dispersion, with u c = for k \\ c, 
whereas for k _L c the plasma frequency is lc^ = (Airne 2 /m^) 1 / 2 . 
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Fig. 1. Spectral functions of La! 9 Sr iCu0 4 at room temperature for ab plane 
(solid curves except (d)) and c axis (dotted curves except (d)). In (d) we display 
uj ab /uj(c) obtained from Eq. (1) with (solid) and without (dotted) correction for the 
phonon-contributions. 



Indeed it is known from optical measurements on cuprates, that with the elec- 
tric held oriented along the planes a plasmon is observed at around 1 eV [1-3] 
whereas above T c there is no zero crossing of Ree c associated with a free-carrier 
plasmon (only optical phonons are observed) [4-7]. In Fig. 1 we demonstrate 
this by plotting the energy loss function (-Iml/e) for both directions as ob- 
tained from our optical experiments. However, the actual situation is more 
complicated in two respects: 

(1) In the normal state the absence of the plasmon along the c direction can not 
be attributed to having m c — > oo, but is rather due to overdamping. This can 
be seen in the following way: With finite damping £ = £00 — u 2 p / u(u -\- i^f ) is the 
(optical) dielectric function at k ~ 0, so that for Au 2 < £qo7 2 the zero crossings 
are along the imaginary frequency axis, i.e. for 2iu = 7 ± {^ 2 — Au 2 / , £ o} 1/ ' 2 • 
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On the other hand, one can use the /-sum rule 



7T 



to obtain u c from the optical data by integrating the imaginary part of the 
dielectric function upto a suitably chosen cutoff frequency x below the on- 
set of interband transitions (in practice such a clean separation is not always 
possible). In Fig. 1(c) we display u p as obtained from Eq. (1) as a function of 
cutoff frequency x along the horizontal axis. The step at 200 cm -1 in the c-axis 
data, as well as a somewhat broader step at 150 cm 1 in the a6-plane data 
are phononic contributions. After correction for these phonon-contributions 
we find that the ratio of electronic plasma frequencies iL> a b/uj c drops from 10 to 
4 as a function of increasing cutoff frequency. The ratios with (solid line) and 
without (dashed curve) correction for the phonon contributions are displayed 
in Fig. 1(d). From LDA band structure calculations we obtained iL> a b/uj c = 5, 
in good agreement with the experimental sum-rule analysis. By fitting the op- 
tical data to a Drude model, we obtained for the optical scattering-rate that 
Tab/Tc > 50, which is consistent with the DC conductivity ratio. 
(2) In the superconducting state there is a partial recovery of a plasmon in 
the c direction, essentially because there exists a finite probability for pair- 
hopping in the superconducting state. This is equivalent to having a small 
but finite critical current perpendicular to the layers, which has a finite pen- 
etration depth associated with it. Using again sum rule arguments, we ex- 
pect that the penetration depth is now given by the dirty limit formula 
2tt\ c = c(3£' 3 cr c (0)//i)~ 1/ ' 2 , which is determined by the 'missing area' in the 
optical conductivity. The zero crossing of Ree occurs now at frequencies below 
the superconducting gap, and corresponds to a collective mode of the super- 
conducting condensate at a frequency uj 2 = c 2 /X^e^. The ratio between this 
frequency and the undamped normal-state plasmon in the ab direction is 



<^ab 



£a6oo^a6arctan (2.3E g /h~{) 
e coo m c 7r/2 



1/2 



The left side of the expression is 50/12000 = 0.009, whereas the righthand side 
becomes 0.011, assuming that E g = 3.5A;bT c and using the other parameters 
as discussed above. Hence if we take into account that m c /m a b < 40 and that 
Tab/Tc > 50, we obtain a consistent description of the anisotropics in cr(0), 
the penetration depth, the screened plasma frequency, and the bare plasma 
frequency obtained from the /-sum rule. 

From the above discussion, we may conclude that the LEG model for the 
anisotropic dielectric properties of the high-T c cuprates is not realistic. The 
major flaw is the fact that taking m c — > oo is a gross oversimplification. The 
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second problem arises due to the strong damping terms. Microscopically the 
anisotropic damping is quite interesting, and has been predicted as a result of 
spin-charge separation [8,9]. To derive the full k- and cu-dependent dielectric 
function including the damping effects is a very difficult theoretical problem. 
The presence of strong anisotropy makes the solution of this problem already 
quite formidable even if the only source of scattering is (elastic) impurity scat- 
tering, which becomes more difficult if the damping has a many-body origin. In 
this paper we present a formula for the k- and u- dependent dielectric function 
without taking into account damping. Although there is no hope that such 
an expression for e gives a complete description of the optical properties or 
electron-energy loss function of high-T c cuprates, the expressions are never- 
theless quite relevant in the discussion of the optical sum rule. Even though 
the optical spectral shape will be very different from that of the undamped 
and unrenormalized plasma poles discussed above, the imaginary part of the 
dielectric function should still obey Eq. (I). 

Let us extend the a layered electron gas model[10,ll], by including a finite hop- 
ping between the planes. The electrons are confined to a stack of 2D sheets 
a distance l c apart, coupled by a finite interlayer hopping parameter t. The 
dispersion relation close to the Fermi-level is 

E k = E F + v F h(k\\ - k F ) + 2t cos(kJ c ) (3) 



We indicate such systems as a coupled layered electron gas (CLEG). The 
plasmon-dispersion was derived for such a model by Grecu[I2] using pertur- 
bation theory. The bare polarization propagator for ksT — > and q — > can 
be expressed as an integral over the Fermi surface 

n°(g» = 2^/ *,f (4) 

87r J J hto — hq ■ v(s) 

s F 



Using straightforward mathematical manipulations, the full expression for 
the k and u dependent polarizability can be easily obtained. To simplify 
the expression we introduce the effective velocity perpendicular to the lay- 
ers v c = 2h~ 1 tl c . In calculating LI° we use the fact that the angular averages 
of terms of the type cos sin 2n+1 are all zero, so that we obtain 

4:TT 3 nv F l c J J u — [ypqw cos <p — v c qj_ sin 0) 



In reciprocal space the expression for the Coulomb interaction between elec- 
trons confined to an infinite stack of ^-layers is v q = ^ire 2 S\q\~ 2 where S 
is the form-factor of the single electron wavefunctions. For ^-layers this is 
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S(ql c ) = S(x) = 1 \x\ 2 sinh(x||)/{cosh(x||) — cos(xj_)}. The dielectric func- 
tion as obtained from the random phase approximation is e = 1 — v q Yl°. In the 
limit q — > we obtain the plasma frequencies w p \\ = e(2kpv p) 1 1 2 (hi c )~ l I 2 and 
u P ± = 2et(2kpl c ) 1 ^ 2 h~ 3 ^ 2 v F 1 ^ 2 for the directions of propagation parallel and 
perpendicular to the planes. In the case of LEG the dielectric function is 



e{q,to) 



1-2S 



LtJ, n 



\q\v F 



qf\V 2 F 



-1/2 



or 



(6) 



On the other hand, from LDA band structure calculations we obtained, that 
for doped La 2 Cu04 the ratio u v i_lu v \\ 0.2, from which we conclude that the 
full plasmon-dispersion formula should be used. The CLEG dielectric function 
using the RPA approximation is 

e(u,q) = 

l2 y,»(!£!l,^) j^o^a] (7) 

where p° is the function 

Zi {' 1/2 

p°(a, b) = — / {(1 - a cos </>) 2 - b 2 }~ cos cjidcj) (8) 
ira J L J 
o 



Let us first consider the plasma dispersion. In Fig. 2 we display this disper- 
sion for the LEG and the CLEG models. The dispersion of the upper branch 
calculated with these parameters is in agreement with electron energy loss ex- 
periments of high-T c cuprates[13,14]. The c-axis plasma frequency is the lower 
bound of a continuum of plasmon states. In the limit q — > both S(ql c ) — > 1 
and p° — y 1, and the optical conductivity for E parallel and perpendicular to 
the planes becomes o rj = Upj/An^j — iu). 

Let us first point out that the above formula for e contains two nontrivial 
terms. In the small-g limit (where the distinction between longitudinal and 
the transverse dielectric functions disappears), the first term describes the 
dielectric function for fields parallel to the planes, whereas the second one en- 
ters the description for electric fields perpendicular to the planes. Although we 
have not been able to give a microscopic justification hereof, the experimental 
results [4,5,15,1] indicate that the many-body effects affect these two terms 
in a very different way: in the cuprates there is a strong anisotropy (of the 
order of hundred), both in u 2 ^ and in the scattering rate ^ r In fact, the c-axis 
conductivity turns out to be practically frequency- independent, whereas the 
in-plane conductivity has a relatively small (but frequency dependent) scatter- 
ing rate. These phenomena cannot be explained from weak elastic or inelastic 
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Fig. 2. Plasmon dispersion of a layered electron gas (left) and that of a coupled lay- 
ered electron gas with = 0.2. The reciprocal-space coordinates are indicated 
as [<7||/ c , qj_l c ] at a number of points. For the Fermi velocity we assumed v F jl c LO\\ = 1. 



scattering of fermionic charge carriers [16,17], but rather require charge car- 
riers with a renormalization factor Z which is very small or zero. Anderson 
and coworkers address all these issues within a single framework of carrier 
confinement to the planes due to the formation of a 2D Luttinger liquid, with 
a de- confinement due to the onset of superconductivity [8,9]. 
As we discussed above, the effect of finite scattering is a quite formidable the- 
oretical problem yet to be solved. However, from inspection of Eq. (7), and 
Fig. 2, we see that the hrst term determines the dispersion of the plasma 
frequency with g||, whereas the dispersion with q± due to the second term is 
relatively small. A simple, and somewhat crude 'shortcut', which is only valid 
for the collective- mode part of the dielectric function (i.e. for frequencies out- 
side the particle hole continuum), is to replace each of the two terms with a 
g-dependent plasma pole, each of which has a different damping term. This 
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suggests that we can make the following approximation in the normal state 
\q\ 2 u[uj + il s {q,LO)\ \q\ 2 u 



This approximation is consistent with our experimental results at k=0. The 
modifications of this function in the superconducting state are quite small, 
and require a very thorough analysis, as in e.g. [18,19]. This dispersion and 
the corresponding loss functions could be measured with e.g. electron energy 
loss spectroscopy or with inelastic X-ray scattering. 
Conclusions 

From an analysis based on the /-sum rule we showed, that the experimental 
values for the anisotropy in plasma frequency along and perpendicular to the 
planes is between 4 and 10. We conclude that the absence of plasmon 
in the normal state of high-T c cuprate superconductors results from a non- 
Drude smearing of the optical conductivity perpendicular to the CuC>2 planes 
over a very wide energy range. This assignment alleviates the discrepancy by 
two orders of magnitude existing in the literature, between the plasma fre- 
quency obtained from optical experiments and that from the single-electron 
interlayer-hopping matrix element. We present an analytical expression for 
the k and u dependent dielectric function for an electron gas, with a Fermi 
surface in the shape of a corrugated cylinder. We propose a phenomenological 
extension of this expression, to take into account the incoherent nature of the 
charge transport in the c direction. 

Acknowledgments We gratefully acknowledge Ir. H. S. Somal for his assis- 
tance with the experiments, and V.H.M. Duijn, N.T. Hien and A. A. Menovsky 
for providing single crystals. This investigation was supported by the Nether- 
lands Foundation for Fundamental Research on Matter (FOM) with finan- 
cial aid from the Nederlandse Organisatie voor Wetenschappelijk Onderzoek 
(NWO). 



References 



[1] J. H. Kim, I. Bozovic, J. S. Harris Jr., W. Y. Lee, C.-B. Eom, T. H. Geballe, 
and E. S. Helfman, Physica C 185-189 (1991) 1019. 

[2] D. van der Marel, J. van Elp, G. Sawatzky, and D. Heitmann, Phys. Rev. B 
37 (1988) 5136. 

[3] S. Tajima, G. D. Gu, S. Miyamoto, A. Odagawa, and N. Koshizuka, Phys. 
Rev. B 48 (1993) 16164. 

[4] Jae Kim, B.J. Feenstra, H.S. Somal, Wen Y. Lee, A. M. Gerrits, A. Wittlin, 
and D. van der Marel, Phys. Rev. B 49 (1994) 13065. 



7 



[5] Jae H. Kim, H.S. Somal, D. van der Marel, A.M. Gerrits, A. Wittlin, V.H.M. 
Duijn, N.T. Hien and A. A. Menovsky, Physica C 247 (1995) 297. 

[6] K. Tamasaku, Y. Nakamura, and S. Uchida, Phys. Rev. Lett. 69 (1992) 1455. 

[7] J. Schiitzman, S. Tajima, S. Miyamoto, and S. Tanaka, Phys. Rev. Lett. 73 
(1994) 174. 

[8] D. G. Clarke, S. P. Strong, and P. W. Anderson, Phys. Rev. Lett. 72 (1994) 
3218. 

[9] S. Chakravarty and P. W. Anderson, Phys. Rev. Lett. 72 (1994) 3859. 
[10] A. Fetter, Ann. Phys. 88 (1974) 1. 

[11] H. Morawitz, I.Bozovic, V.Z. Kresin, G. Rietveld, and D. van der Marel, Z. 
Phys. B 90 (1993) 277. 

[12] D. Grecu, Phys. Rev. B 8 (1973) 1958. 

[13] N. Nuecker, H. Romberg, S. Nakai, B. Scheerer, J. Fink, Y. F. Yan, and Z. X. 
Zhao, Phys. Rev. B 39 (1989) 12379. 

[14] N. Nuecker, U. Eckern, J. Fink, and P. Mueller, Phys. Rev. BAA (1991) 7155. 

[15] D. van der Marel, M. Bauer, E. H. Brandt, H.-U. Habermeier, D. Heitmann, 
W. Konig, and A. Wittlin, Phys. Rev. B 43 (1991) 8606. 

[16] A. G. Rojo and K. Levin, Phys. Rev. B 48 (1993) 16861. 

[17] M. J. Graf, D. Rainer, and J. A. Sauls, Phys. Rev. B 47 (1993) 12098. 

[18] D. van der Marel, Phys. Rev. B 51 (1995) 1147. 

[19] Y. Zha, K. Levine, and D. Z. Liu, Phys. Rev. B 51 6602 (1995). 



8 



